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We model the dynamics of attractively interacting ultracold bosonic atoms in a quasi-one-dimensional wave¬ 
guide with additional harmonic trapping. Initially, we prepare the system in its ground state and then shift the 
zero of the harmonic trap and switch on an additional narrow scattering potential near the center of the trap. After 
colliding with the barrier twice, we propose to measure the number of atoms opposite to the initial condition. 
Quantum-enhanced interferometry with quantum bright solitons allows us to predict detection of an offset of 
the scattering potential with considerably increased precision as compared to single-particle experiments. In a 
future experimental realization this might lead to measurement of weak forces caused, for example, by small 
horizontal gradients in the gravitational potential — with a resolution of several micrometers given essentially 
by the size of the solitons. Our numerical simulations are based on the rigorously proved effective potential 
approach developed in [Phys. Rev. Lett. 102, 010403 (2009) and Phys. Rev. Lett. 103, 210402 (2009)]. We 
choose our parameters such that the prerequisite of the proof (that the solitons cannot break apart, for energetic 
reasons) is always fulfilled, thus exploring a parameter regime inaccessible to the mean-field description via the 
Gross-Pitaevskii equation due to Schrodinger-cat states occurring in the many-particle quantum dynamics. 
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1. INTRODUCTION 

Matter-wave bright solitons can be used for interferometry 
by profiting from their mean-field description based on the 
Gross-Pitaevskii equation (GPE) — primarily by the fact that 
attractive interactions prevent the wave-packets from spread¬ 
ing [1-5]. Matter-wave bright solitons' are investigated ex¬ 
perimentally both for attractive interaction [3, 7-15] And — 
in the presence of an optical lattice — also for repulsive inter¬ 
action [16]. Recent experimental results for bright solitary 
waves in attractively interacting Bose-Einstein condensates 
include both a first step to interferometric applications [3] and 
collisions of two bright solitons [12] as well as quantum re¬ 
flection off a narrow attractive potential barrier [13]. 

While current experiments can successfully be modeled 
by using the mean-field (GPE) approach, beyond-mean-field 
properties of quantum bright solitons are a focus of ongo¬ 
ing theoretical investigations [17-28]. In the current paper, 
we explore beyond mean-field many-particle quantum prop¬ 
erties of matter-wave bright solitons [29, 30] (cf. [27]) for 
interferometric purposes. While for mean-field based inter¬ 
ferometry, quantum fluctuations in, for example, the center of 
mass velocity can endanger the interferometric scheme [31], 
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' Our solitons are not solitons in the strict sense but rather solitary waves. 
However, these solitary waves can behave very similar to solitons [6]. 


in the present paper both interactions and quantum fluctua¬ 
tions are the tools to generate quantum-enhanced interferom¬ 
etry. While spin-squeezed states [32-37] are one possibility to 
achieve [38] quantum-enhanced interferometry, Schrodinger- 
cat states^ are another [47-49] (cf. [50]). An experimental re¬ 
alization of an atom interferometer can be found in Ref. [51]; 
interferometry using non-entangled states leading to enhanced 
precision can be found in Ref [52]. 

We use the rigorously proved [53] effective potential ap¬ 
proach developed independently of each other in Refs. [21, 
22]. The effective potential approach is valid in the regime of 
low kinetic energies where the soliton is energetically forbid¬ 
den to classically break apart, thus corresponding to a mean- 
field regime where reflection at a barrier leads to steps in the 
transmission coefficient [29, 54, 55].^ In this energy regime 
both mesoscopic Schrodinger-cat states [21, 29] and, for two- 
component Bose condensates, mesoscopic Bell states [58] 
have been predicted theoretically to occur in the many-particle 
quantum dynamics of quantum matter-wave bright solitons. 
Eor a single-species, attractively interacting Bose condensate 
and repulsive barrier, the effective potential approach [21, 22] 
is a powerful tool to model the many-particle quantum dy¬ 
namics both for solitons of the order of N - 100 atoms [21] 
but also for much smaller particle numbers of N - 2 [59]. 


^ Schrodinger-cat states are discussed in Refs. [21, 23, 30, 39-46] and refer¬ 
ences therein. 

^ In the regime of higher kinetic energy, where mean-field bright solitons 
do break apart when hitting a barrier, scattering has been investigated, for 
example, in Refs. [4, 56, 57] and references therein. 
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The paper is organized as follows: Section II describes 
how quantum bright solitons in one-dimensional wave guides 
can be modeled; followed by the effective potential approach 
(Sec. Ill) we use to describe the many-particle quantum dy¬ 
namics in the regime where matter-wave bright solitons are 
energetically forbidden to break apart. In Sec. IV we in¬ 
troduce the signature we propose to use for future experi¬ 
ments: the probability to find all particles on the side op¬ 
posite to the initial condition after the soliton hits the bar¬ 
rier for the second time. In Sec. V we present our results: 
the quantum-enhancement achieved by using quantum bright 
solitons rather than single particles (or non-interacting Bose- 
Einstein condensates). The paper ends with conclusions and 
outlook in Sec. VI. 


II. MODELING QUANTUM BRIGHT SOLITONS 

In order to model attractively interacting atoms (gin < 0) in 
one dimension, the integrable Lieb-Liniger-(McGuire) Hamil¬ 
tonian [60, 61] 




2m dxj 


N-l N 

ZZ glDS(Xj - Xn) 

j=\ n=j+\ 


( 1 ) 


is available, where xj denotes the position of particle j of mass 
m. The (attractive) interaction 


giD = 

< 0 


( 2 ) 


is proportional to the s-wave scattering length a and the per¬ 
pendicular angular trapping-frequency, cuj_ [62]. Including 
the center-of-mass momentum K, the (internal) Lieb-Liniger 
ground state [19, 61] reads (cf. [19]) 


'P(xi, X 2 ,..., Xiv) oc exp 


m\g\ 


2n^ 
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the center-of-mass coordinate is given by 


(3) 


1 ^ 
N 




The ground-state energy is given by 

r,2 Ml AA 


Eo(N,giD) 


rngj^NiN^ - 1) 
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(4) 


(5) 


Even more relevant for what we propose here, there is an en¬ 
ergy gap between the internal ground state and the first excited 
state. 


£gap = Eo(N - 1) - Eo(N) 

_ mg^io^iN- 1) 

“ 8^2 ■ 


( 6 ) 

( 7 ) 


As long as the total kinetic energy of the matter-wave bright 
soliton lies well within this energy gap. 


Ekin ^ Ep 


( 8 ) 


the validity of the effective potential approach of Refs. [21,22] 
can be proved rigorously [53]. This parameter-regime coin¬ 
cides with a complete break-down of the validity of the Gross- 
Pitaevskii equation"^ for describing scattering matter-wave 
solitons off, for example, a narrow barrier [29, 54, 55, 58]. 

If the center-of-mass wave function is a delta function and 
the particle number is A » 1, then the single-particle density 
can be shown [19, 64] to be equivalent to the mean-field result 
based on the Gross-Pitaevskii equation [63] (cf. footnote 4) 


g(x) = 


N 


(9) 


4^iv(cosh[x/(2^iv)]] 

normalized to the total number of particles A; the soliton 
length is given by 


= 


m\g\D\ (N - 1) 

Adding a longitudinal harmonic trapping potential. 


( 10 ) 


N 

y(Xi,X 2 , ...,Xiv) = y VoC^c;), (11) 

1=1 

Vo(x) = ^mo?{x + Xof, (12) 

does not change the physics (beyond breaking integrability) 
- as long as [65] the soliton length is small compared to the 
harmonic oscillator length 



The plane waves in the center-of-mass coordinate in Eq. (3) 
become harmonic oscillator eigenfunctions in this limit [66]. 
After preparing the system in its ground state, at f = 0 the 
single-particle potential Vo is changed quasi-instantaneously 
to 

Viix) - ^ma/x^ + vi6{x - As), (14) 

thus moving the center of the longitudinal trap by Aq and 
adding a very narrow scattering potential (modeled by a delta 
function) which is shifted by a distance of As from the trap 
center. As the proof of the validity of the effective potential 
approach introduced in the following section uses repulsive 
potentials [21, 53] we choose vi > 0, where vi quantifies the 
strength of the scattering potential. 


The mean-field approach via the Gross-Pitaevskii equation (GPE) [63] 
d fP' 

in—tip = - — —ip + {N- i)giD|p|2p. 

is a useful tool to explain ongoing state-of-the art experiments with matter- 
wave bright solitons like the experiments of Refs. [12, 13]. However, the 
mean-field approach fails in the quantum regime [29, 58] that makes the 
quantum-enhanced interferometry proposed in the cuiTent paper possible. 
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III. THE RIGOROUSLY PROVED [53] EFFECTIVE 
POTENTIAL APPROACH 


The effective potential method was developed indepen¬ 
dently by K. Sacha, C. A. Muller, D. Delande, and J. Za- 
krzewski [22] and C. Weiss and Y. Castin [21]. The rigorous 
mathematical proof by Y. Castin published a couple of years 
later [53] includes strict error bounds demonstrating that this 
is a highly reliable many-particle method in the energy regime 
where bright solitons cannot break apart for energetic reasons 
[Eq. (8)]. 

While this restricts the applicability to kinetic energies that 
lie within the energy gap between the internal ground state 
(all N particles in one soliton) and the internal first excited 
state {N - 1 particles in one soliton), the effective potential 
is a high-end many-particle method that replaces the many- 
particle Schrodinger equation by an effective equation for the 
center-of-mass wave function in a non- perturbative way [53]. 
In cases where the center-of-mass wave function is localized 
at two places separated by more than the soliton length, this 
automatically corresponds to highly entangled Schrodinger- 
cat states. 

In general, the effective potential is the convolution be¬ 
tween the scattering potential and the soliton [21, 22, 53]. For 
the relevant single-particle potential (14) the result is partic¬ 
ularly simple; the harmonic trap simply is replaced by a har¬ 
monic trap with identical trapping frequency but for a parti¬ 
cle of mass niN', the delta function is replaced by the mean- 
field soliton profile (9) multiplied by the strength of the delta- 
function scattering potential vi (cf. [53]). As both potentials 
are a function of the center-of-mass coordinate, this leads to 
a further huge advantage of the effective potential approach: 
by introducing the effective potential, we have replaced the 
many-particle Schrodinger equation by an effective single¬ 
particle Schrodinger equation for the center of mass X with 


in—ip(x,t) 

at 




^{X,t) 


2mN dX^ 

+ VMip(x,t), 
1 


Veff(Y) ^-niNcJx^ 


Ekin ^ Egap, 

viN 


(15) 


4^^(cosh[(Y-Ys)/(2^w)]] 


2 ■ 


(16) 


The proof of the effective-potential approach requires that the 
kinetic energy is smaller than the energy gap (7) [21, 53]. 

For sufficiently small ratio of soliton length to harmonic 
oscillator length, the effective potential can even be approxi¬ 
mated by [29] 


Veff(Y) ^ ]^mN(j?X^ + vxN5{X - As), (17) 

which could be treated analytically [67]. In the following sec¬ 
tion, our approach will be to combine an approximate analyt¬ 
ical treatment of the scattering of the soliton off the potential 
with full numerics for which the center-of-mass wave function 
of the ground state soliton is prepared in the harmonic oscilla¬ 
tor ground state. The trap is then shifted such that the soliton 


moves towards the new minimum — in which then a narrow 
scattering potential is switched on. 

While the model displayed in Eqs. (15)-(17) looks like a 
single-particle approach, these equations for a particle of mass 
Nm describe many-particle quantum dynamics in the param¬ 
eter regime in which quantum bright solitons cannot break 
apart [21, 53] (cf. [22]). Any solution of Eq. (15) which 
spreads over more than the soliton width can thus be iden¬ 
tified as being a mesoscopic quantum superposition [21] of 
the Schrodinger-cat type relevant for quantum-enhanced in¬ 
terferometry [47]. In the following, we will focus on 50:50 
splitting leading to the particularly useful [47] high-fidelity 
Schrodinger cat states. To identify more general quantum su¬ 
perpositions see, e.g.. Ref. [68]. 

To summarize, thanks to the rigorous proof of Ref. [53], the 
effective potential approach is a highly reliable many-particle 
method for quantum bright solitons. Its level of reliability 
in the regime of low kinetic energies (in which the proof of 
Ref. [53] is valid) is at least as high as high-end many-particle 
methods used, for example, in [23, 26] (cf. Refs. [69, 70]). 
The occurrence of highly entangled Schrodinger-cat states rel¬ 
evant for quantum-enhanced interferometry [47, 49] is thus 
guaranteed — under conditions that are within the reach of 
ongoing state-of-the-art experiments [11-14]. 


IV. EXPERIMENTALLY ACCESSIBLE SIGNATURE 


We propose to prepare the ground state of N attractively 
interacting atoms in a one-dimensional wave-guide with addi¬ 
tional week harmonic trapping. After quasi-instantaneously 
shifting the harmonic potential by several soliton lengths, 
we switch on a very narrow repulsive scattering potential in 
the center of the trap and observe the quantum dynamics 
for 50:50 splitting after the first collision for two collisions 
(cf [29, 30, 45]). In Fig. 1 we compare what happens to a 
single atom {N - 1) to a A = 100 quantum bright soliton; In 
the limit that the effective potential behaves similar to a delta- 
function, the length-scales are a factor of 10 smaller for the 
A = 100 soliton^ compared to A = 1. In passing we note 
that the signature demonstrated in Fig. 1 b (and also in 3 c) 
is not accessible on the mean-field (GPE) level as it involves 
intermediate Schrodinger-cat states® as the soliton is not en¬ 
ergetically allowed to break into smaller parts, thus what is 


^ For ^Li and N » 100, the set of parameters used is the slightly modified 
parameter set (we doubled the strength of the radial trapping frequency) of 
Ref. [21] for the i-wave scattering length a = -1.72x 10“® m, cj± = IttxIx 
4800Hz. In addition we chose w = IttxSHz. The initial distance from the 
trap center is chosen to be Xq = 30 ^m; this can easily he reduced for bright 
solitons (but not for single particles or weakly interacting Bose-Einstein 
condensates). For the soliton it corresponds to iikiii/|fto(iV)| ^ 0.83, within 
the regime of validity of the effective potential approach [cf. Eq.(17)]. 

® In the energy regime in which the many-particle quantum mechanics pre¬ 
dicts Schrodinger-cat states, the GPE shows strong jumps [29, 54, 55] in 
the transmission behavior. When colliding with the scattering potential 
twice, the solitons are either reflected or transmitted in both cases. Thus, 
the GPE-soliton ends on the same side as the initial condition (cf. [56]) and 
on the opposite side to the many-paificle quantum prediction. 
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FIG. 1. (Color online) As the experimentally accessible signature we 
suggest to use the probability to find all particles on the side oppo¬ 
site to the initial condition after scattering twice [shown is the two- 
dimensional projection of the modulus squared of the center-of-mass 
wave function as a function of both position and time for (a) A = 1 
and (b) a A = 100 soliton for the parameters of footnote 5]. After 
being prepared in the ground state of a harmonic trap, at r = 0 the 
trap is shifted and a narrow repulsive scattering potential switched 
on at the new center of the trap. After the scattering potential is hit 
for the first time, Schrbdinger-cat states occur, after the second colli¬ 
sion all particles are on the side opposite to the initial condition (cf. 
[29]). The fact that Schrbdinger-cat states occur is a consequence 
of the fact that slow bright solitons cannot energetically break apart 
(cf. [21], footnote 5); the interaction remains switched on the entire 
time. 


shown in Fig. 1 b for times t ^ 100 ms corresponds to a quan¬ 
tum superposition of all particles being either on the left or 
right of the scattering potential. The main difference between 
the two panels of Fig. 1 is that the upper panel is only for one 
particle. While one might choose to call this a Schrodinger 
cat [39], if describing a non-interacting Bose-Einstein con¬ 
densate with N particles in a rather classical product state, the 
following section shows that the state displayed in the upper 
panel lacks the many-particle entanglement properties neces¬ 


sary for quantum-enhanced interferometry for which the state 
displayed in the lower panel is particularly useful. 

In order for the signature displayed in Fig. 1 to work, the 
time-scale for decoherence events via, for example, single- or 
three- particle losses has to be large compared to the total time 
of the experiment from the first time the soliton hits the barrier 
to being detected at the side opposite to its initial condition. 
Reference [21] shows a possible parameter set for ^Li, using 
conservative estimates for the loss rates (cf. [71]), thus justify¬ 
ing that using trapping frequencies in the direction of motion 
of the order of w = 27r x 10 Hz is indeed justified. 

Before we show how the signature displayed in Fig. 1 can 
be used for future experiments with quantum-enhanced mea¬ 
surements, we would like to stress again that the fact that we 
have Schrodinger cat states relies on strong attractive interac¬ 
tion preventing the soliton from breaking apart and thus allow¬ 
ing the rigorously proved [53] effective potential approach of 
Refs. [21, 22] Eqs. (15)-(17) to be used. While often interac¬ 
tion is only present to generate entanglement, for the signature 
displayed in Eig. 1 to work it is furthermore essential to leave 
the interaction switched on the entire time — otherwise the 
feature that all particles end up opposite to their initial condi¬ 
tion would not work. 


V. QUANTUM-ENHANCED MEASUREMENT: 

BEYOND-CLASSICAL PRECISION 

A. Overview of Sec. V 

In this section, we use the mathematically rigorous effec¬ 
tive potential approach (Sec. Ill) to show that we can obtain 
beyond-classical precision in interferometric measurements 
with quantum-bright solitons. The fact that we are in the 
regime of validity of the effective potential approach guaran¬ 
tees both the occurrence of Schrodinger cat states and ensures 
that this is not accessible to non-linear physics (for the GPE 
in the regime of low kinetic energies, the soliton would either 
pass the barrier twice or be reflected twice, thus in both cases 
ending back in the original position rather than opposite to the 
initial position as shown for the quantum case in Eig. 1. 

In Sec. V B we show that the signature to find all parti¬ 
cles opposite to the initial condition oscillates as a function 
of the distance of barrier from the center of the trap (moved 
by a horizontal force) position [Eq. (20)] both by using an 
intuitive analytical approximation which is backed up by sys¬ 
tematic numerical simulations. In Sec. V C we summarize the 
role Schrodinger cats play in our scheme and stress the role of 
measuring atoms in Bose-Einstein condensates via scattering 
or absorbing photons. Section V D shows how this is related 
to the precision; as for the quantum-enhanced interferometry 
with Schrodinger-cat states described in Refs. [47, 48] gain 
a factor of N in the precision with which the phase can be 
detected. Repeating the single-particle experiment N times 
more often than the experiment with A-particle Schrodinger- 
cat states yields an overall win of a factor of 1/ V(V for the 
Schrodinger cats [47, Sec. “Quantum-Enhanced Parameter 
Estimation”]. In Sec. VE we show how well horizontal dif- 
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ferences in the gravitational field could be measured. 


B. Dependence of the signature on the distance of the harrier 
from the trap center 

As a practical example showing how much this improves 
the precision, we perform a thought measurement of a hor¬ 
izontal gradient in the gravitational potential, modeled as a 
linear potential added to Eq. (15), 

l^gravity(A^) ~ ( 18 ) 

which would lead to a shift of the zero of the harmonic trap¬ 
ping potential and thus to a non-vanishing distance 




gravity 


(19) 


between trap center and scattering potential. In order to es¬ 
timate the precision of our interferometric measurement, we 
need to determine the position of a maximum by repeating the 
measurements n times. The interference pattern as a function 
of Xg — namely the transmission coefficient after scattering at 
the potential twice — is given by the approximate analytical 
formula [29] 


T = 
Q = 


- [1 H- COS(4fIXs)] 
NmcoXo 

n 


( 20 ) 

( 21 ) 


In order to derive the last line, we used energy conservation to 
derive the center-of-mass momentum fiK via 


{TiKf 


1 


= -NmaP'Xi, and thus 


2Nm 
K = 


2 

NmcoXo 

n 


( 22 ) 

(23) 


together with the text-book [72] result that scattering a particle 
of mass Nm with a plane wave with momentum HK from a 
delta function potential 


y(A) = —5(X) (24) 

Nm 

leads to fifty-fifty splitting if 

K^Q.. (25) 

While the true effective potential will be broader than the 
delta-function used here, this only affects the amplitude, not 
the spatial period of the oscillations (see Fig. 2 a, cf [29, 
Fig. 3]). For the same parameters as the numerical data de¬ 
picted in Fig. 2 b. Fig. 2 c shows the analytical formula (20) 
confirming that this formula predicts the width of the interfer¬ 
ence patterns correctly. As predicted by Eq. (21), the spatial 
period increases for decreasing Xq (see Fig. 2 b). Thus the 
delta-function approximation is sufficient for the purpose of 
estimating the precision of our interferometric setup. 




FIG. 2. (Color online) (a) Our signature that everything is on the side 
opposite to the initial condition after colliding with the barrier twice, 
denoted by T, as a function of distance between trap minimum and 
scattering potential for different widths of the 1 / cosh^-potential, in¬ 
troduced in Eq. (16), Xg = 30 pm and N = 100 particles. The value 
^ = 0.42 pm corresponds to the parameters from footnote 5. The 
limit ^ > 0 corresponds to a delta function potential, (b) Same for 

varying initial displacements Xg and f = 0.42 pm kept fixed, (c) Dis¬ 
plays the analytical result (20) for the same parameters as panel (b). 
Thus panel (c) confirms that the analytical approximation (20) qual¬ 
itatively correctly capture the physics displayed in the many-particle 
quantum dynamics of panel (b). 


C. Role of Schrodinger-cat states & measuring particles in a 

BEC 


The signature we use for our interferometry scheme re¬ 
lies on a very fundamental level on Schrodinger-cat states: 
after two collisions with a barrier, neither the numerical re¬ 
sults for the transmission after two collisions with a batTier 
depicted in Fig. 2 nor the analytic result of Eq. (20) could 
have been obtained with a mean-field description via the GPE. 
In the regime of low kinetics energies in which the rigor¬ 
ously derived effective potential approach is valid (see Sec. Ill, 
[21, 22, 53]), the GPE predicts unphysical jumps in the trans- 
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mission/reflection behavior [29, 54, 55]; solitons would be ei¬ 
ther reflected twice or transmitted twice, in both cases end¬ 
ing up at the same side of the scattering potential for the two 
collisions investigated here, corresponding to the mean-field 
(GPE) prediction 


7’gpe - 0. 


(26) 


Rewriting Eq. (20) to show the A^-dependence explicitly. 


T 


1 + cos N 


^4ma»2£'o2fs 

n 


(27) 


yields a reduction of the width of the interference patterns by 
a factor of N. This is identical to the general case discussed 
in Ref [47, Sec. “Quantum-Enhanced Parameter Estimation”] 
when replacing a single-particle quantum superposition by an 
A^-particle Schrodinger cat, the first step towards quantum- 
enhanced interferometry. 

Observations of atomic Bose-Einstein condensates employ 
scattering or absorption of laser light [73]. Multiple scattering 
let alone absorption of the same photon are not an issue that 
is discussed in such cases [73], a fact that is further stressed 
by the small atom numbers 100 which are the focus of 
our current paper. Thus, measuring a BEC of N noninteract¬ 
ing atoms is equivalent to measuring N single atoms indepen¬ 
dently. In both cases, the central limit theorem thus increases 
the precision by a factor of 1 / ^/N. 

Note that if measuring the Schrodinger cat state, an indi¬ 
vidual measurement of T is either 0 or N. While having N 
atoms in a bright soliton can simplify detecting it, there is no 
additional correction corresponding to the factor 1 / Viv. 

To determine the precision, in Sec. VD repeat the above 
thought experiment n times (comparing an non-interacting 
BEC with the bright soliton Schrodinger-cat state). 


D. Precision 


In order to answer how precisely we can determine the po¬ 
sition of the scattering potential with respect to the trap center, 
we have to consider a couple of separate points. All are based 
on the distance between minima/maxima, as determined by 
Eq. (20). Half the distance between two minima is given by 


AXs 


n 

4Q 

1 nfi 
N 2majA'o 


(28) 

(29) 


1. The aim is to move the scattering potential back to the 
center when it was moved by a small distance from 
there; thus our device will have to remember by how 
much the potential has moved. 


2. Repeating the experiment n times for various distances 
will give us a part of the cosine of Eq. (20) centered 
around the maximum. 


3. Half the distance between two minima gives a first or¬ 
der of magnitude of the precision with which we can 
measure the position; the larger the particle number in 
our soliton the higher the precision will be. (In an ac¬ 
tual experiment, the full width at half maximum could 
be used equivalently). 

4. If measured in terms of the precision that can be read 
off the distance of the interference fringes, a quantum- 
enhanced measurement is better by a factor of 1/A. 
This is the typical scaling for quantum-enhanced inter¬ 
ferometers with Schrodinger-cat states as compared to 
single particles [48, Fig. 1]. This difference of a factor 
of N is visible by comparing Fig. 3 a and Fig. 3 c. 

5. In our case, the increased precision described in point 4 
above is related to the de Broglie wavelength of the par¬ 
ticles; smaller wavelengths allowing to measure smaller 
distances. Thus, this is independent of precise details of 
the initial state. 

6. Repeating the quantum-enhanced measurement n times 
with Schrodinger-cat states of N atoms in order to de¬ 
termine the position of the maximum involves n x N 
atoms. 

7. In order to compare the outcome with the single-particle 
case (cf. [48]) thus repeat the same interferometric se¬ 
quence of measuring the position of the maximum with 
single particles nxN times. Compared to the quantum- 
enhanced case, this increases the number of measure¬ 
ments by a factor of N, allowing us to determine the 
position of the maximum more precisely. As the errors 
involved in repeating such experiments N times more 
often scale as 1/ VA [47, Sec. “Quantum-Enhanced Pa¬ 
rameter Estimation”], 


non-interacting particle /on\ 
0-^N non-interacting particles ^ ^ 

yN 

we thus have 

non-interacting particles 

OA;y_particle soliton ^ 

Viv 

SXn=1 non-interacting particle 
oc - 

A 

Thus, to summarize the precision of the single-particle case 
is limited by the Heisenberg-uncertainty relation [the distance 
measured (29) is of the order of the de Broglie wavelength; 
independent of the precise initial state, the momentum is of 
the order of matXo]. The minimal distance measurable can 
be improved by statistical analysis of repeating the experi¬ 
ment A times, giving an improvement of a factor of 1 / Va for 
the smallest distance that can be measured [Eq. (30)] for re¬ 
peated experiments with single particles [47, Sec. “Quantum- 
Enhanced Parameter Estimation”]. A precision which scales 
faster with A as displayed in Eq. (32) thus seems, at first 
glance, to violate the Heisenberg uncertainty relation. 


(31) 

(32) 
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However, many-particle entangled states like Schrodinger- 
cat states are known to show such behavior which is some¬ 
times described as “beating the standard quantum limit” [47]. 
As the de Broglie wavelengths are a factor of N smaller than 
for the single-particle case, it should not be suprising that 
we can measure distances with the precision of Eq. (32). 
The overall gain of a factor of 1 / ViV is of course consistent 
with the Heisenberg-uncertainty relation if applied to our N- 
particle Schrodinger-cat state (cf. [47, 48]). 

While the analytical calculations are based on the approx¬ 
imate result given in Eq. (20), the oscillation periods are 
confirmed by doing full numerics (Eigs. 2 and 3): In panel 
Eig. 3 a, the full numerical result shows that it is straightfor¬ 
ward to identify the position of the maximum. Combining 
the central-limit theorem with the analytical approximation 
using Gaussians demonstrates that repeating the experiment 
100 times will lead to an increased precision for the position 
of the maximum of a factor of 10 (Eig. 3 b). While the nu¬ 
merics using the delta-function agrees with the analytic pre¬ 
dictions of Eq. (20), for the 1 / cosh^ potential corresponding 
to the parameters of footnote 5 the amplitude is smaller. 


E. Measuring horizontal differences in the gravitational 
potential 

Thus, our thought experiment can measure differences in 
the gravitational acceleration of at least 

lAggravityl < (33) 

which yields a difference in acceleration of 

,. ,1 

|Aggrav.,yl<^^. (34) 

This yields a difference in the acceleration potential of 

1 ntoh 

22fo|Aggravityl < --, (35) 

A m 

which takes into account the fact that we effectively average 
Aggravity over a distance IXq. Thus, within the definition (35), 
the precision is independent of our choice of Xq. The strength 
of quantum bright solitons lies in particular in the fact that 
one can get much closer to the scattering potential than for 
non-interacting particles, thus exploring shorter length-scales. 

Eor ^Li, N - 100 (cf. [21], footnote 5) and to - In x 5 Hz, 
if our thought experiment involves measuring n solitons, it 
can detect differences smaller than (0.009/ ■\/n) fimmls^. By 
replacing ^Li by we could gain another order of mag¬ 
nitude in precision. While at first glance it might appear to 
be tempting to propose to both further reduce the trapping 
frequency and further increase the particle number, the va¬ 
lidity of Eq. (35) relies on the occurrence of intermediate 
Schrodinger-cat states and thus will work better for not too 
large particle numbers and not too long time-scales. 

Depending on the magnitude of the potential differences 
we intend to measure on micrometer scale (gravity was just 



Xsifim] 




- 0.15 0.0 0.15 

Xs[nm\ 

FIG. 3. (Color online) Our signature that everything is on the side 
opposite to the initial condition after colliding with the barrier twice 
as a function of distance between trap minimum and scattering po¬ 
tential [see. Fig. 1 also for parameters used] for (a) a single particle: 
numerical data [magenta (dark gray) dashed curve] and analytical 
formula (20) (black solid curve) and (b) repeating the experiment 
100 times with single particles [magenta (dark gray) dashed curve] 
compared to the single-particle result (20) (black solid curve) and re¬ 
placing this result by a combination of Gaussians [green (light gray) 
dash-dotted curve], (c) Shows quantum-enhancement by a factor of 
1/VA for a quantum bright soliton of A = 100. Shown are for¬ 
mula (20) (black solid curve), numerical data using a delta-function 
barrier [green (light gray) dash-dotted curves] and the 1 / cosh^ po¬ 
tential corresponding to the parameters of footnote 5 [magenta (dark 
gray) dashed curve]. 


an example), the precision might actually be too large to eas¬ 
ily determine the distance of the scattering potential from the 
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trap center. To work out a clever way to start with reduced 
precision (with non-interacting atoms and/or higher trapping 
frequencies) in order to roughly determine the position and to 
continuously increase the precision could still be an engineer¬ 
ing challenge. Nevertheless, our results clearly indicate the 
feasibility of such a device. 


VI. CONCLUSION AND OUTLOOK 

Similar to the case of the spin-squeezed states realized in 
the experiment of Ref. [38] demonstrating interferometric pre¬ 
cision surpassing what is possible classically, Schrodinger- 
cat states can also be used for quantum-enhanced measure¬ 
ments [47, 48] (cf [49]). Quantum bright solitons generated 
from attractively interacting Bose gases are an ideal system to 
produce such Schrodinger-cat states [21]; state-of-the art ex¬ 
perimental setups with ultra-cold atoms [11-14] could reach 
such a regime in the near future. Ultracold atoms can be ap¬ 
plied to questions not easily accessible to photon experiments 
(cf. [74]). 

In the present paper, we propose to use small, harmonically 
trapped quantum matter-wave bright solitons for quantum- 
enhanced measurements. The signature we propose to mea¬ 
sure — the probability to find all particles at the side opposite 
to the initial condition — is not accessible on the mean-field 
(GPE) level as it involves intermediate Schrodinger-cat states 
(see footnote 6). 

Our numerical simulations are done using the effective po¬ 
tential approach of Refs [21, 22] by taking care that the pre¬ 
requisite of the proof [53] is fulfilled: The solitons are en¬ 
ergetically forbidden to break apart. Thus, if the center-of- 
mass wave function is split 50:50 when the soliton hits the 
barrier for the first time, we are certain to have Schrodinger 


cat states. For the second step to work (that after recombining 
at the barrier the soliton ends at the side opposite to the initial 
condition, it is furthermore essential to maintain the quantum 
superposition [29]. While experimentally observing the mo¬ 
tion of quantum bright solitons in the presence of decoherence 
via atom losses [75] would also be an interesting experiment, 
the quantum-enhanced interferometry will only work if deco¬ 
herence via particle losses happens on time-scales longer than 
the duration of the experiment (cf. [21]). 

When shifting the position of the scattering potential away 
from the center, our signature oscillates. The spatial oscil¬ 
lation period scales oc 1 / -s/n with the number of particles, 
thus providing the “usual” increase in precision by a factor of 
1 / ViV quantum-enhanced measurements can provide. Here, 
it occurs when the position of the maximum is determined via 
repeated measurements. 

A potential experimental realization of such a setup might 
be useful to measure small horizontal gradients in the gravita¬ 
tional force on micrometer scales in an approach complemen¬ 
tary to existing experiments combining Bose-Einstein con¬ 
densates and gravity [76]. 

The data presented in this paper are available online [77]. 
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